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Abstract: We analyze the phases of the Sakai-Sugimoto model at finite temperature and 
baryon chemical potential. Baryonic matter is represented either by 4-branes in the 8-branes 
or by strings stretched from the 8-branes to the horizon. We find the explicit configurations 
and use them to determine the phase diagram and equation of state of the model. The 4-brane 
configuration (nuclear matter) is always preferred to the string configuration (quark matter), 
and the latter is also unstable to density fluctuations. In the deconfined phase the phase 
diagram has three regions corresponding to the vacuum, quark-gluon plasma, and nuclear 
matter, with a first-order and a second-order phase transition separating the phases. We find 
that for a large baryon number density, and at low temperatures, the dominant phase has 
broken chiral symmetry. This is in qualitative agreement with studies of QCD at high density. 
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1. Introduction 

QCD at finite baryon density has a rich phase structure (for reviews see Q §). Naively 
one would expect that at high density, like at high temperature, QCD is in a deconfined 
chiral-symmetric quark-gluon plasma phase. It turns out, however, that new phases appear 
at high density, in which both the chiral symmetry and the gauge symmetry are broken |2[ |j . 
In real QCD with N c = 3 and three light flavors of quarks the dominant phase at high density 
is a color-flavor-locking (CFL) phase Q. With two light flavors of quarks the dominant phase 
is a color-superconductor. At large N c it is believed that these gauge-symmetry breaking 
phases are suppressed, and the dominant phase at high density is a chiral density wave ||, [|, 
in which the chiral symmetry (only) is broken non-uniformly. It appears therefore that QCD 
(both with N c = 3 and at large N c ) at low temperature and high density always has broken 
chiral symmetry. 

These results rely on perturbative calculations in QCD, and analogous models such as the 
Nambu-Jona-Lasinio (NJL) model, at finite density near the Fermi surface, and are therefore 



- 1 - 



limited to values of the chemical potential for which a s {jj) -C 1. At present, lattice QCD 
techniques are unable to deal with a (large) baryon chemical potential (for a review see [j?]]). 

At large N c , gauge/gravity duality is an alternative approach to gauge theory at strong 
coupling (8|. Several recent models have incorporated flavors using probe branes in back- 
grounds dual to large N c Yang-Mills theories with various amounts of supersymmetry M. 



The Sakai-Sugimoto model in particular is quite similar to QCD at large N c [1C]. This model 
builds on Witten's model for pure Yang-Mills theory in four dimensions, which uses 4-branes 
wrapped on a Scherk-Schwarz circle [11], and adds Nf probe 8-branes and Nf probe anti-8- 



branes transverse to the circle. These provide massless chiral fermions (left-handed from the 
8-branes, right-handed from the anti-8-branes) in the fundamental representation of both the 
gauge group U(N C ), and the flavor group U(Nj)l x U(Nf)n. 

One of the most compelling features of this model is that it describes spontaneous chiral- 
symmetry breaking in a simple geometrical way. Since in the near-horizon limit the circle 
vanishes at a finite radial coordinate, the 8-branes and anti-8-branes are smoothly connected 
into a U-shaped configuration with an asymptotic separation L at infinity. The actual embed- 
ding of the 8-branes is determined by solving the DBI equations of motion with this boundary 
condition. 

The model also exhibits many other properties similar to QCD [12, [HJ. In particular it 



has an interesting phase structure at finite temperature [14]. At low temperature the model 
is essentially the same as at zero temperature, i.e. it describes a confining gauge theory with 
broken chiral symmetry. At high temperature the model deconfines and chiral symmetry 
is restored, which is described geometrically by the separation of the 8-branes and anti-8- 
branes. For sufficiently small L there is also an intermediate range of temperatures at which 
the model is deconfined but chiral symmetry remains broken. In the deconfined phase both 
the connected U-configuration and the separated parallel configuration are possible. The 
dominant configuration, and therefore phase, is determined by comparing their actions. 

The baryonic U(l)y symmetry corresponds in models with fundamental matter to the 
diagonal U(l) gauge symmetry of the probe branes. Baryon number is therefore described by 
electric charge, and baryon number density is related to the electric field, or more precisely to 
the electric displacement field, of the diagonal U(l). Correspondingly, the baryon chemical 
potential is described by the value of the gauge potential at infinity ^4o(°°)- Finite baryon 
density in the Sakai-Sugimoto model has been studied in [15, [l6|, [l?], |i8j . However, only part 



of the parameter space has been explored so far. Other models with finite baryon density 
have been studied in |J, ||, |l], ||. 

In this paper we explore the full parameter space of the Sakai-Sugimoto model at finite 
temperature and finite uniform baryon number density, in both the confined and deconfined 
phases. Other than temperature and baryon number density (or chemical potential), this 
model has an additional parameter not present in QCD, namely the asymptotic 8-brane-anti- 
8-brane separation L. We will assume that the value of L is such that the intermediate phase 
of deconfinement with chiral symmetry breaking exists, in other words that the deconfine- 
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ment temperature is (much) smaller than the chiral-symmetry restoration temperature. The 
confined phase is, of course, of great interest, but the deconfined phase exhibits a much richer 
phase diagram. The deconfined phase is also qualitatively the same as the non-local NJL 



model [23, 24], and we expect a similar phase diagram in that case. 

There are two types of objects which carry baryon charge. The baryons themselves 
correspond to 4-branes wrapped on the 5 4 part of the background. Due to the RR flux each 



4-brane comes with iV c strings attached 25]. The other end of the strings is attached to 
the 8-branes, which is how the baryons get their flavor. However, in the nonsupersymmetric 
4-brane background this configuration is not static. The strings pull the wrapped 4-branes up 



toward the 8-branes [26]. When they reach the 8-branes, the 4-branes can also be described as 
instantons in the world-volume theory of the 8-branes. In the deconfined phase baryon charge 
can also be carried by strings which stretch from the 8-branes all the way to the horizon. 1 
This describes a possible phase in which baryon charge is carried by free quarks. 

In both cases we will consider a uniform distribution in M 3 of baryon charge, so the 8- 
brane worldvolume theory reduces to a one-dimensional problem in the radial coordinate with 
a source term. As the 4-branes and strings exert a force on the 8-branes, their embedding at 
finite baryon number density will have a cusp. 

We will make the following approximations. First, we will assume that the wrapped 4- 
branes are pointlike in the transverse coordinates and uniformly distributed in M 3 . The precise 
description would be in terms of instantons in the 8-brane worldvolume theory. However, in- 
stanton solutions in DBI are not known. Analysis in the Yang-Mills (plus Chern-Simons) 
approximation shows that an instanton has a finite size on the order of the string length |l3j . 
Using pointlike instantons, or equivalently pointlike 4-branes, is therefore a good approxima- 
tion. We will also neglect any direct interactions between the 4-branes themselves or between 
the strings. 

Our results are summarized in the phase diagram in figure |]. At low temperatures the 
theory confines, and there is a second order phase transition at finite [i to a phase of nuclear 
matter. At high temperature the theory deconfines and chiral symmetry is restored. At 
intermediate temperatures chiral symmetry is broken for all ji (as in QCD). We find a second 
order phase transition also in the intermediate temperature range between the vacuum and 
nuclear matter phases. This is similar to QCD, but in QCD it is a first-order transition due 
to the attractive interaction between the baryons (which we have neglected). 

In section 2 we begin by describing the possible 8-brane configurations corresponding 
to the Sakai-Sugimoto model at finite baryon number density. In section 3 we discuss the 
thermodynamics of the gauge theory which are implied by these configurations and derive 
the full phase diagram in the grand canonical ensemble. We consider both the confined and 
deconfined phases, and both 4-branes and strings as sources of baryon charge in the deconfined 
phase. It will turn out that 4-branes are always preferred to strings, and that the stringy 
"quark matter" phase is actually unstable to density fluctuations. We conclude and offer 



1 We can think of these strings as ending on massless, wrapped 4-branes at the horizon. 
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Figure 1: The phases of holographic QCD at finite temperature and baryon chemical potential. A 
particular deconfmcmcnt temperature (0.025) was chosen for illustration purpose only. 

suggestions for future work in section 4. 

2. Finite density brane configurations 

The basic brane configuration consists of Nf 8-branes and Nf anti-8-branes in the near 
horizon background of N c 4-branes wrapped on a Scherk-Schwarz circle with N c S> Nf. At 
zero temperature the background is capped and the circle is topologically trivial, so the 8- 
branes and anti-8-branes connect into a U-shaped configuration. The dual gauge theory is 
confining, and chiral symmetry is broken. At finite temperature this continues to be the 
only possible configuration until one reaches a critical temperature, at which the dominant 
background switches to the black hole, and both U-shaped 8-branes and separated parallel 
8-branes and anti-8-branes are allowed. At first the U-shaped configuration dominates, so 
chiral symmetry remains broken even though the gauge theory deconfines. Chiral symmetry 
restoration occurs at a second critical temperature, which for L small enough, is above the 
first critical temperature (otherwise they are equal). At this temperature the separated 8- 
brane-anti-8-brane configuration begins to dominate. 

The baryon number current is related holographically to the diagonal £7(1) part of the 
8-brane gauge field. To study finite baryon number density configurations we therefore need 
to include this gauge field in the 8-brane action. The first place it enters is in the DBI action: 



Sds = — A*8 j d 9 X e * Tr \f-det(g M N + ^ol'Tmn) 



(2.1) 



where T is the U{Nf) field strength 



T = dA + iA A A . 



(2.2) 
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We decompose the U(Nf) gauge field into an SU(Nf) part and a U(l) part as follows 2 

A = A + —^A. (2.3) 

The U(l) gauge field will also appear in the CS action, which will be important below. Let's 
study the effect of turning on this gauge field on the brane configuration in the different 
phases. 

2.1 confined phase 

In the confined phase the background (at finite temperature) is given by 

3 3 

ds 2 = 5 {{dXif + (dX) 2 + f{U)dXl) + * (j*L + (2.4) 

fu\ 3/4 



r) 

F 4= ^f\ (2.6) 

where ~ X§ + (3, X 4 ~ X 4 + /3 4 , and 

t/^ „ /47r\ 2 i? 3 



/(CO = 1 - %F > = (y ) ^ , i? 3 = TtgsNM'? 12 ■ 
It is convenient to express everything in terms of dimensionless quanitities, so we define 



(2.7) 



U X 4 Xf A 27TQ'i . , 

u = t: > x 4 = tt > r = - ?r > a = — ; • (2-8) 

i? R R ^/2N]R v ; 



The 8- brane action with the U(l) gauge field is then given by 3 

S m =M J duu 4 f(u)(x' A {u)f + ± (-L _ (a' ( U )) 2 ) 



(2.9) 



where we have defined the overall normalization as 

(2 . 10) 

9s 

where ^4 is the volume of a unit S 4 , and V3 is the volume of space (M 3 ). Note that the action 
scales as NfN c . 



2 We are using the convention Tr T a Tb — ^8 ab . 

3 This is the action for just the 8-branes, i.e. for 1/2 of the full configuration. The lower limit of the integral 
is the lowest radial position of the 8-brane configuration, and the upper limit is infinity. 
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As will become clear in the next section it is convenient to also define the Legendre- 
transformed action 



Sd8 = Sd8+N / dud(u)a (u 



where d(u) is the electric displacement field defined by 

1 5S D8 ua 
d ( u ) = ~ 77 



M Sd' Q (u) 



/(u)(xi(u))2 + «-s (-^-(fi^n)) 



(2.11) 



(2.12) 



This gives 



Sds = A/" / c?n u 



/(n)(4(n)) 2 + 



n 3 /(n) 



1 + 



(d(«)) s 



n" 



(2.13) 



The equations of motion for x^{u) and d(u) can be integrated once yielding two constants: 

d(u) = d 



(x' 4 (u)f 



f(u)(u* 



fid 2 ) 



1 



-1 



(2.14) 



J(uoM+u 3 d 2 ) 

where no is defined as the position where x 4 (n) diverges. 

For d = the solution is a U-shaped 8-brane in the (x4,n) plane, with no as its lowest 
radial position (figure |2[) However, a non-trivial electric displacement d requires a source 
at u = n c , which is possibly different from no, which will change the 8-brane configuration. 
This is essentially a one-dimensional electrostatics problem in the coordinate n, except that 
the 8-brane covers the region [n c ,oo] twice. Each part carries an electric displacement d. 

The only possible sources for d in the confined phase are instantons, or equivalently 4- 
branes wrapped on the S 4 , in the 8-branes. For a uniform d we need a uniform distribution 
in M 3 of 4-branes. The source term comes from the 8-brane CS action: 



Scs 



Ms 
6 



R^xR+xS 4 



C 3 Tr {iTra'Tf 



24vr 2 



R*xR + 



UJ 5 (A). 



The relevant term is the one that couples the U(l) to the SU(Nf): 



24^ 2 J RixR+ y/W } 

We will assume a uniform distribution of 4-branes in 

1 



3 - 9 
yl TrF 2 . 



at n = u c '. 



— -TrF 2 = n4<5(n — n c )d 3 xdn. 
87r^ 



(2.15) 



(2.16) 



(2.17) 



where is the (dimensionless) density of 4-branes wrapped on S 4 . The equation of motion 
for the U(l) gauge field then gives 



d'{u) 



2Tia'R 2 M 



714 5(u — U c ) 



(2.18) 
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Figure 2: The 8-brane configuration with d = and d ^ in the confined phase. 



and therefore 4 



rt4 



(2.19) 



The instanton distribution ( |2.17| ) also sources the equation of motion for X4(u) and will 
therefore deform the shape of the 8-brane. Physically, the 4-branes pull down on the 8-branes. 
Since the 4-brane distribution has a finite energy density per unit 7-volume (the S 4 they wrap 
plus the M 3 ), it will form a cusp in the 8-brane (like a bead on a string). Away from the 
cusp the 8-brane will follow two opposite pieces of a U-shaped solution, which are truncated 
at some radial position u c above uq (figure ||). 

The value of u c can determined by the zero- force condition in the (x4,u) plane. The 
proper tension of the 8-brane is given by varying the Legendre-transformed action ( |2.13 ) 
(which is the same as the Hamiltonian once we substitute in for the solution of £4(11)) with 
repsect to the proper distance along the 8-brane. The result is 



2\ 1/2 



(2.20) 



The force due to the 4-branes is given simply by varying their action with respect to their 
position u c , again taking care to vary with respect to the proper distance. The 4-brane action 
is 



Sd4 = ra4 ^^ 4 J dn 4 dT e * y 7 detg M N 



-Mu c d , 



and the force is therefore 



dS 



D4 



1 



TD4 



du c Jg. 



till: 



\Mdul' A ^J{^). 



(2.21) 



(2.22) 



u=u c 



4 This is the density of 4-branes on one half of the 8-brane configuration. The total 4-brane density is twice 
this much. 
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Figure 3: The position of the cusp (and 4-brane) in the 8-brane as a function of the electric displace- 
ment d for I = 1 in the confined phase. We present this as a log plot to show both the initial decrease, 
as well as the limiting value at large d. 



The condition for equilibrium is then 

f D 8Cos6 = f m , (2.23) 
where is the proper angle of the 8-brane at u c , 



cos 9 



f( U0 ) (tig + tigd 2 ) 1 2 



(2.24) 



An elegant alternative derivation of this result is given in the Appendix. 

We want to solve this for u c , while holding fixed the asymptotic separation of the 8-branes 
and anti-8-branes, which is given by 

1 = 2 dux' 4 (u). (2.25) 

Ju c 



This can be done numerically by varying uq and d, computing u c and I using ( 2. 23] ) and 



( |2.25| ), and then tabulating (d,u c ) for a given value of I. The result is presented in figure 
||, where for definiteness we have set 1 = 1. Note that for small values of d the cusp comes 
down as d increases, but beyond a certain value of d it goes up. Initially the 4-branes pull 
the 8-branes down, but eventually the 8-branes win this "tug-of-war" . We will see the same 
behavior in the deconfined phase. The initial downward motion of the cusp indicates that 
the chiral condensate of the gauge theory decreases initially as the density increases, and the 
eventual upward motion indicates that the chiral condensate eventually increases with d. As 
we will soon see this has an important implication for the gauge theory thermodynamics at 
high density. 

2.2 deconfined phase 

The background describing the deconfined phase is given by exchanging the roles of X4 and 



t: 



ds 2 = J {f{u)dr 2 + (dx) 2 + dx\) + «-§ (j^- + u 2 <ml\ , (2.26) 
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Figure 4: Possible 8-brane configurations with d = and d ^ in the deconfined phase. 



with the same dilaton and RR 4-form as before, and where 



4vr\ 2 i? 2 /4tt\ 



(2.27) 



where t = R/ (3 T = RT is the dimensionless temperature. The 8-brane action is now given by 
S m =Af f duu 4 [f(u)(x' 4 (u)) 2 + u- 3 {l-(a' (u)) 2 )]^ , (2.28) 



and the Legendre-transformed action is 

S DS =M j duu A [f{u){x' A {u)f + u- 3 ] 
where d(u) is now given by 

d(u) 



1 + 



(d(u)) 



2lh 



u;' 



(2.29) 



ud' (u) 



[/(u)(xi(u))=» + «-s(l-(a&(u))a)]s 



(2.30) 



As in the confined phase, the equation of motion for d(u) implies that it is a constant d(u) = d. 
On the other hand, for £4(14) there are two types of possible configurations (figure |j). The 
first corresponds to separated parallel 8-branes and anti-8-branes with 



x' 4 (u) = . 

and the second to a connected configuration with 



(2.31) 



1 



u 3 f(u) 



f(u)(u 8 + u 3 d 2 ) 



(2.32) 



Actually, as we will soon see there are in fact two connected solutions, but only one is 
(classically) stable. 



1 

. 8 
. 6 




Figure 5: The position of the cusp in the 4-brane cusp configuration as a function of the electric 
displacement for I = 1 in the deconfined phase. 



2.2.1 4-brane sources 

The parallel configuration can have a uniform electric displacement d without sources. In the 
connected configuration, however, we need a source, as in the confined phase. One possible 
source is again 4-branes inside the 8-branes. The 4-brane action is now 

1 



S D 4= gAAneVT^Jd, (2-33) 



and the force they exert is given by 



Sm = \Nd (/(«.) + ^) uT = \Md^^uT , (2.34) 
where in the last equality we used the form of f(u) in (2.27). The 8-brane force is computed 



as before, but with the metric of the deconfined phase. This gives 

, / ,72x1/2 

f m =Mufl A ^J{^)[l + -^ , (2.35) 

and the same angle as before ( [2.24| ). The solution of the zero- force condition for a representa- 
tive temperature is shown in figure ||. The qualitative behavior is the same as in the confined 
phase: initially the cusp comes down as d increases, but eventually it goes up and approaches 
a fixed value. 

In the deconfined phase there are actually two connected solutions in general. This can 
be seen by looking at I as a function of the cusp position u c at fixed d and t (figure ||). There 
are two values of u c for a given I below some l max . At I = l max the two solutions coincide, 
and above l ma x there is no connected solution. This behavior is also true for d = 0. So 
in fact there are three solutions in all when I < l ma x'- the parallel configuration, a "short" 
cusp configuration (or "short" U-connguration when d = 0), and a "long" cusp configuration 
(or "long" U-connguration when d = 0). When I > l ma x only the parallel configuration is 
a solution. This picture is very reasonable from the following point of view. Imagine that 
we have two classically stable solutions in some theory with a potential. The two solutions 
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Figure 6: The asymptotic brane-anti-brane separation I as a function of the cusp position u c for a 
fixed d = 0.5 and t = 0.1. 




Figure 7: A schematic of the three possible solutions: (a) I < l max (b) I = l max (c) I > l ma x- 

correspond to two local minima of the potential. But this necessarily implies that there should 
be a third solution, corresponding to the local maximum between the two minima (figure |7]a). 
This solution should be unstable. Now imagine that one local minimum is lower than the 
other, and that the second local minimum approaches the local maximum as we vary some 
parameter. When they coincide we get a point of inflection (figure 0b). As we continue to 
vary the same parameter both solutions cease to exist, leaving only the lower minimum (figure 
0c). This is precisely what happens for the 8-brane embedding. The parallel and short cusp 
configurations are the stable solutions. The long cusp configuration must therefore correspond 
to the unstable solution. We leave it as a future excercise to exhibit the required negative 
mode. Note that this picture necessarily implies that we don't have to worry about the 
cusp solution disappearing, since in the region of parameter space near this point the parallel 
solution always dominates. 

2.2.2 string sources 

The other possible sources of electric displacement in this phase are strings which stretch 
from the 8-branes to the horizon at uj\ 5 We can determine the precise relation between the 

5 Instead of ending at the horizon, the strings may also end on 4-branes which wrap the S 4 and are located 
below the 8-branes. This is Witten's description of baryons Each 4-brane has iV c strings attached to it. 
However this configuration is not a solution of the equations of motion. There is a net force that pulls the 
4-brane to larger u Eventually the 4-brane reaches the 8-brane and turns into an instanton. 
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Figure 8: The position of the cusp in the stringy cusp configuration as a function of the electric 
displacement for 1 = 1. 



density of strings n s and the electric displacement d by looking at the i?-neld dependence of 
the supergravity, 8-brane, and string actions: 

Ssug[B) = / d w x^/^U)e- 2 *\dB\ 2 (2.36) 

S D8 [B] =M J duu 4 \j{u){x' 4 {u)f + n- 3 - u~ 3 (B 0u + a^(u)) 2 ] 5 (2.37) 
Sfi[B) = J drdu (y-det g M N ~ B 0u ) . (2.38) 

Varying with respect to Bq u and integrating over an 8-sphere surrounding the endpoint of 
the strings in the 8-branes we find that 

27Ta'R 2 M 

n s = ^7 d. (2.39) 



Note that this is consistent with what we found for 4-branes in ( |2.19| ), since each 4-brane 
(away from the 8-branes) has N c strings attached. 

Evaluating the string action for the deconfined background gives 

S F1 = N{u c - u T )d . (2.40) 

As in the 4-brane case, we have assumed a uniform distribution of strings in M 3 x S 4 , so the 
point on the 8-brane where they end will again be a cusp. The force downward applied by 
the strings is given by 



SSfi 1 

JFl 



5u c Jgl, 



AAdn^vTK). (2.41) 



The solution to the zero-force condition with the strings is shown in figure ||. The behavior 
is different from the 4-brane case. The position of the cusp comes down monotonically with 
increasing d. 

It turns out, however, that the stringy cusp configuration is always subdominant to the 
4-brane cusp configuration. We can see this by comparing their actions. The total action will 



- 12 - 





Figure 9: Comparing the actions (relative to the parallel configuration) of the string-sourced (red) 
and 4-brane-sourced (blue) cusp configurations. The 4-brane case wins at all temperatures and all d. 



have a contribution from the 8-branes Sds given by ( 2.2S| ), where the integral is taken from 
u c to infinity, and from either the 4-brane or string sources. The integrals are divergent, but 
we can regularize them by subtracting the action of the 8-branes in the parallel configuration 
in both cases. The results are shown in figure ||. The action of the cusp configuration sourced 
by 4-branes is smaller than that of the configuration sourced by strings at all temperatures 
and for all values of d. 

In the next section we will show that in fact this configuration is unstable to fluctuations in 
d. This is similar to the instability found in [19]. We will comment on a possible interpretation 
of this instability in the conclusions. 



3. Thermodynamics with finite chemical potential 

We now turn to the gauge theory implications of the configurations we found. Our main 
goal is to understand the phase diagram of the gauge theory at finite temperature t and 
finite baryon chemical potential \x. Note that this model has an additional parameter not 
present in QCD corresponding to the asymptotic 8-brane-anti-8-brane separation I. We will 
generally fix I = 1. We have also considered other values of I (smaller and larger) and found 
no qualitative change in the results. 

3.1 baryon chemical potential 

The grand canonical potential is obtained by evaluating the 8-brane action ( |2.9[ ) or ( |2.28D on 
the solution. For convenience we will normalize the potential by dividing out the normaliza- 
tion constant J\f, 

fi(t, /■*) = S D8[t, Xa{u), d (u)] solution ■ (3.1) 

Note however that the potential, as well as all other thermodynamic quantities associated to 
the matter, scale as NfN c . The (dimensionless) baryon chemical potential /i is identified with 
the asymptotic value of the U(l) gauge potential in the solution 



H = ao(oo) 



(3.2) 



With our normalizations the baryon number density is given by 6 



n b 



dfjL 



d. 



(3.3) 



We will therefore use d to denote also the density. 

For computational purposes it is more convenient to express \i in terms of d using the 
canonical ensemble. The free energy is defined as 



F(t,d) = n(t,n) + lid, 
and the chemical potential is given by 



dF(t, d) 



dd 



(3.4) 



(3.5) 



The free energy is thus related to the Legendre-transformed 8-brane action on the solution. 
In the cusp configuration the total free energy includes also the contribution of the source 
4-branes or strings, evaluated at the position of the cusp: 



F(t,d) = — -. (jS D& [t,X4(u),d(ll)] solution + S source (t, d,u c 



(3.6) 



The dependence on d comes from three places: the explicit dependence of Sds and S source on 
d, the dependence of the solution for x'^ on d, and the dependence on d of u c . Including all 
of these gives 

solution 



(3.7) 



The second term vanishes since 8Sd%/8x'± (u) is constant by the equation of motion for X4 and 
the integral of 8x' 4 /8d, at fixed u c , gives 81/ dd which vanishes since I is fixed. The third and 
fourth terms cancel by the zero- force condition at the cusp (see Appendix), leaving 




f°° 1 

L + A7 



dS„ 



dd 



(3.8) 



t.l.Uc 



where o,q(u) is related to d by inverting the relation ( 2.12j ) or ( 2.3Q ). The identification of the 
chemical potential with the value of the gauge potential at infinity (3.2) therefore reflects a 
particular choice of gauge, in which clq{u c ) is identified with the mass of the source. In the 
parallel configuration the source term vanishes, and the lower limit of the integral is at the 
horizon u = ut- In this case the gauge choice (|3.2| ) gives o,q{ut) = 0, which is consistent with 
the fact that the source becomes massless at the horizon. 



6 The true baryon number density is given by (p.l 
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3.2 confined phase 

In the confined phase only a connected 8-brane configuration is possible. However for a given 
//, that is at fixed do(oo), there are two connected solutions, a U-configuration with d = 
and a 4-brane sourced cusp configuration with d ^ 0. The former corresponds to the QCD 
vacuum, and the latter to a phase of nuclear matter. In the vacuum phase do is constant, 
and the electric displacement d vanishes. Therefore O does not depend on fj, in this phase. 
In the nuclear matter phase ao(u) is sourced by 4-branes, and the chemical potential is given 
by (|3.8|) , which in the confined phase yields 

f°° d 1 

H = / du — ^^^^^^^^^^^=^^^^^^^^^^^ H — u c . (3-9) 



lie 



y/f(u) ( U 5 + (i 2 ) _(M) 3 /(uo) {ul + d*) 



3 



Note that both u c and uq depend on d and I. There is no temperature dependence in the 
confined phase. As we are working in the grand canonical ensemble, this represents an implicit 
expression for d(fi). We see that there is a critical value for the chemical potential 2u c /3 for 
which d = 0. Below this value there is no cusp solution and therefore no nuclear matter 
phase. This is precisely the onset chemical potential ^onset- For (j, > \i nset both the vacuum 
and the nuclear matter phases exist, and we must compare their grand canonical potentials 
to determine which phase is preferred. These quantities are actually divergent at u — > oo, 
but the difference is finite and is given by 

u 5/2 u 5/2 



u 5/2 



(3.10) 



where uo n and uq v are the (different) values of no in the nuclear and vacuum phases, respec- 
tively. The results, for the representative value ukk = 0.5 (and I = 1), are shown in figure 
|ic| . Since AQ < for all u > u onse4 , the nuclear matter phase is preferred. Figure [l0| also 
shows the baryon number density as a function of the chemical potential for ukk = 0.5. (The 
value of Honset grows with ukk)- Near \x = [i nset the density is continuous and behaves as 
d ~ (u — fJ-onset) 1 - This is therefore a second-order phase transition with a critical exponent 
of 1. This is a reasonable result given our approximation of ignoring the baryon interactions. 



3.3 deconfined phase 

In the deconfined phase the situation becomes more interesting since there are more possible 
configurations at a given value of u. In addition to the U-shaped and 4-brane-cusp config- 
urations, there are the parallel configuration, with vanishing or non-vanishing density, and 



- 15 - 




Figure 10: The difference in the grand canonical potential fl between the nuclear and vacuum state, 
and the resulting baryon number density, as a function of /i in the confined phase. 



the string-cusp configuration. The parallel configuration corresponds to a phase in which 
chiral symmetry is restored. At finite density this is the quark-gluon plasma (QGP). The 



vacuum parallel configuration is irrelevant, since it is clear from ( 2.28| ) that a non-trivial d is 
always preferred in the parallel configuration. The string-cusp configuration features strings 
stretched between the 8-branes and the horizon. We will refer to the corresponding phase in 
the gauge theory as quark matter. However, as we shall soon see this configuration is actually 
unstable, at least for a uniform distribution of baryon charge. That leaves three phases to 
compare: the vacuum (U-configuration) , nuclear matter (4-brane-cusp) and the QGP (finite 
density parallel configuration). We therefore expect in general a phase diagram in the 
plane with three phase regions. 

3.3.1 unstable quark matter 

Let us first show that the quark matter phase (string-cusp configuration) is thermodynam- 
ically unstable. Evaluating the chemical potential ( |3.8| ) in the deconfined phase with string 
sources ( |2,40D yields 



H= I du Vf(uj d + ( Ue _ nT ). 



Uc 



y/f(u) (^ + d 2)-(^) 3 /M (tig + d 2 ) 



Figure 11 shows a plot of d vs. fi for this configuration. It is apparent that 

|<0, (3.12) 

(or equivalently dfi/dd < in the canonical ensemble) and therefore that the string-cusp 
configuration is thermodynamically unstable to density fluctuations. A similar instability 



was found in the D3-D7 model at finite density [19|. 



3.3.2 phase diagram 

We begin by comparing the vacuum phase to the QGP phase. The vacuum phase is described 
by the U-configuration, and the quark-gluon plasma (QGP) phase is described by the parallel 
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Figure 11: Baryon number density vs. chemical potential in the string-cusp configuration in the 
deconfined phase. 



configuration. The grand canonical potential of the vacuum can be read off from (2.28) and 
( g3g ) with d = 0: 



^vac V 



M=r*> / /2 ^ • (3.i3) 

^f(u)-^f(u ) 



The potential of the QGP is given by ( jjQp and ( p3o| ) with a^(u) = 0: 



n 9flP (^) = ^ dn^===. (3.14) 

The density d is a function of /z which is obtained from Q3.8| ) without sources, and ( p.30| ). 
This gives 

f°° d 

/ du , (3.15) 

Ju T vu° + d z 

which can be inverted numerically to get d(^). Both potentials are divergent at u — ► oo but 
the difference is finite: 

Af^ = Qqgp — Q V ac • (3.16) 

Figure ^ shows AOi(/i) for a few representative temperatures (and 1 = 1). The transition 
between the two phases occurs when AQi = 0. We find a line of transitions between these 
two phases in the (t, fi) plane shown in figure 12. This result was obtained previously in [16]. 



We now turn to the comparison of the vacuum phase with the nuclear matter phase. We 
did this already in the confined phase, and found a second-order phase transition at some 
A 1 = Honset- Since the Sakai-Sugimoto model exhibits chiral-symmetry breaking also in the 
deconfined phase, it is reasonable to expect that nuclear matter should form also in this case. 
The potential of the nuclear phase (with 4-branes) is given by 



fWM = / du VA ' = , (3.17) 
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Figure 12: Grand canonical potential and phase diagram for the vacuum vs. QGP phases. 



o . 002 



-0 . 002 
-0 . 004 
-0 . 006 
-0 . 008 





. 1 


0T2\ 


. 3 




. 1 



0.2 0.3 



Figure 13: Grand canonical potential and baryon number density in the nuclear matter phase relative 
to the vacuum phase. 



where d is again given implicitly in terms of /j, using (|3.8|) with 4-brane sources: 



Ju ° V/(«) (« 5 + d 2 ) - 3 /(V) K + rf 2 ) d 

We are now interested in the difference between the potentials of the nuclear phase and the 
vacuum phase, 

A$^2 = S^nuc — Qvac • (3.19) 

Figure |l^ shows Af^A*) for a representative temperature. The behavior is qualitatively 
the same at all temperatures: A^2 is negative for all \i for which the nuclear phase exists. 
Figure [l^ also shows the density as a function of fj, at the same temperature. As in the 
confined phase, the critical exponent is 1 to within our numerical accuracy, so the transition 



is second order. By varying the temperature we obtain the phase diagram in figure 14. The 
behavior agrees qualitatively with what is expected in QCD: ^onset decreases slightly as the 
temperature increases. 

The final part of the phase diagram comes from comparing the nuclear and QGP phases: 

ASI3 = fi nuc — Qqgp • (3.20) 
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Figure 14: Phase diagram for vacuum and nuclear matter phases. 
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Figure 15: Grand canonical potential (for t — 0.1,0.12,0.15) and phase diagram for the nuclear vs. 
QGP phases. 



Here we find an interesting temperature dependence. Figure 15 shows AQ^(pt) for three 
representative temperatures. At low temperature the nuclear matter phase wins for all /x. 
Then there is a temperature range for which the system undergoes two transitions as /j, is 
increased, first from nuclear matter to QGP, and then back to nuclear matter. The resulting 
phase diagram is shown in figure 15. The physical source of the dip in the phase diagram is 
the dip that occurs in the position of the cusp u c as a function of the density d (figure ||). 
There is a similar dip in the phase diagram of QCD (see for example Q). At high temperature 
(not shown) the QGP phase is preferred for all \i. 

Finally, combining the three separate phase diagrams gives the complete phase diagram 
shown in figure 1. At low temperature and chemical potential the vacuum phase dominates, 
at low temperature and high chemical potential the nuclear phase dominates, and at high 
temperature chiral symmetry is restored and the quark-gluon plasma phase dominates. 



3.4 entropy and equation of state 

Phases of thermodynamic systems are also characterized by the their equation of state and 
entropy. Let us briefly discuss these for the different phases we have encountered. 

The pressure as a function of the density p(d,t) is essentially given by — Q(fj,(d),t). We 
find that at low temperature the behavior in the confined and deconfined phases is very 
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Figure 16: Density vs. chemical potential for large densities. 



1.5 



0.5 



0.05 0.1 .15 T 



Figure 17: Entropy vs. temperature in the deconfmed phase. 

similar. 7 At small densities d ~ (fx — fJ- on set) and therefore 8 

p(d) ~ (m - Honset) 2 ~ d 2 . (3-21) 

At large densities d ~ /x 5 / 2 (figure UJ]) and thus 

p(d) ~ pJl 2 ~ d 7 / 5 . (3.22) 

It is interesting that although we have not specified that the baryons are fermions (indeed 
there seem to be both fermionic and bosonic components), the results for n(d) mimic a 
behavior expected for fermions. This is due to the response of the DBI action to the electric 
field. 

The entropy as a function of the temperature s(t) is computed from the free-energy 
F(t, d). The interesting case is the deconfmed phase, since there is no temperature dependence 
in the confined phase. At low temperature, where chiral symmetry is broken, we find (for 
both small and large densities) 

s(t) ~ t 5 , (3.23) 
and at the high temperature, where chiral symmetry is restored, we find 

s(t) ~ t 6 . (3.24) 



7 The behavior at high temperature, i.e. in the QGP phase, was essentially worked out in p4| . 
3 For a free fermi gas p(d) ~ d 5//3 . 
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4. Conclusions 

In this paper we have analyzed the different phases of the Sakai-Sugimoto model at finite 
temperature and baryon chemical potential and determined the phase diagram. In many 
respects our phase diagram is similar to that of QCD. In both cases chiral symmetry is 
broken at low temperature and restored at high temperature, at all values of the chemical 
potential. The dip in the phase diagram suggests that the chiral condensate initially decreases 
with n and then increases. This is similar to the behavior in QCD. It would be interesting to 
study this directly using the holographic description of the chiral condensate in terms of the 
tachyon p7|| . 

The finite density phase is described by a gas of "baryonic matter" 4-branes wrapped 
on the S 4 inside the 8-branes. At low temperatures this phase always dominates over the 
quark-gluon plasma phase. The other possible description of baryon matter in terms of strings 
("quark matter") turns out to be subdominant and unstable. If we ignore these facts and use 
strings instead of 4-branes to describe baryonic matter the phase diagram would change and 
chiral symmetry would be restored at high density. 

We also find a phase transition between the vacuum and nuclear matter phases. In QCD 
this is a first-order transition, but in our case it is second-order. We believe that that the 
difference is a result of neglecting the interactions between the 4-branes. 

Another difference is that QCD at high density is expected to be in a CFL phase, in 
which both the chiral symmetry and the gauge symmetry are broken. However, at large iV c 
QCD is expected to be dominated by a non-uniform chiral symmetry breaking phase with 
unbroken gauge symmetry. We did not explore this possibility in this paper. However the 
result that the "quark matter" phase, with strings stretched to the horizon, was unstable to 
density fluctuations, suggests that there may exist a stable non-unifrom phase. It would be 
interesting to see if it is similar to the chiral density wave in large N c QCD. 

A. Zero-force condition from the action 

The force balance condition for the cusp configurations can alternately be obtained directly 
by varying the total action with respect to the cusp position u c . The total action is given by 



We want to vary the action with respect to u c , while keeping the physical varaibles I, t and d 
fixed. To do this we need to vary x' 4 (and therefore uq) accordingly. This gives 







However since I is given by 




(A.3) 
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we get 

du d A 



-x' A {u c ) + / 

J U r . 



0. (A.4) 



d,tl 



Furthermore, the equation of motion sets 5Sd&/$x' 4 to a constant independent of u. Requiring 
the total action to be stationary with respect to the varaition of u c one gets 

t i \ ii \$Sd8 9S source . 

L W -x 4 (« c )^_. (A.5) 

Substituting in the expressions for L and S source then reproduces the force balance conditions 
in the various cases (confined, deconfined, 4-branes, strings). 
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